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We present measurements of the branching fraction, the polarization parameters 
and C-P-violating asymmetries in — > D*^D*~ decays using a 140 fb^^ data sam- 
ple collected at the T(45) resonance with the Belle detector at the KEKB energy- 
asymmetric e+e- collider. We obtain B{B^ D*+D*-) = [0.81 ± 0.08(stat) ± 
O.ll(syst)] X 10-3, i?_L = 0.19 ± 0.08(stat) ± O.Ol(syst), Rq = 0.57 ± 0.08(stat) ± 
0.02(syst), S = -0.75 ± 0.56(stat) ± 0.12(syst) and A = -0.26 ± 0.26(stat) ± 
0.06(syst). Consistency with Standard Model expectations is also discussed. 

Key words: B decay, CP violation, sin2^i 
PACS: ll.30.Er, 12.15.Ff, 13.25.Hw 



1 Introduction 

In the Standard Model (SM), CP violation arises from an irreducible complex 
phase, the Kobayashi-Maskawa (KM) phase [1], in the weak- interaction quark- 
mixing matrix. In particular, the SM predicts CP asymmetries in the time- 
dependent rates for 5° and decays to a common CP eigenstate fcp [2]. 
Recent measurements of the CP- violation parameter sin20i by the Belle [3,4] 
and BaBar [5] collaborations established CP violation in — > J/ipK^ and 
related decay modes [6], which are governed by the b — > ccs transition, at a 
level consistent with KM expectations. Here 0i is one of the three interior 
angles of the Unitarity Triangle [3,4]. 

Despite this success, many tests remain before it can be concluded that the 
KM phase is the only source of CP violation. The S° 0*^0*" decay, 
which is dominated by the h — > ccd transition, provides an additional test of 
the SM. Within the SM, measurements of CP violation in this mode should 
yield the sin 20i value to a good approximation if the contribution from the 
penguin diagram is neglected. The correction from the penguin diagram is 
expected to be small [7]. Thus, a significant deviation in the time-dependent 
CP asymmetry in these modes from what is observed in 6 — > ccs decays would 
be evidence for a new CP- violating phase. 

In the decay chain T(45') B^B^ — > /cp/tag, where one of the B mesons 
decays at time tcp to a final state fcp and the other decays at time itag to 
a final state /tag that distinguishes between P° and the decay rate has a 
time dependence given by [2] 



g-|At|/rso 



{ 



l + q Ssm{AmdAt) + Acos{AmdAt) 



} 



(1) 



4tbo 
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Here S and A are CP-violation parameters, tbo is the 5° lifetime, Ama is 
the mass difference between the two mass eigcnstates, At = tcp ~ ttag, 
and q = +1 (—1) when the tagging B meson is a B^ (B^)- The parameter S 
corresponds to the mixing-induced CP violation and is related to sin 20i, while 
A represents direct CP violation that normally arises from the interference 
between tree and penguin diagrams. 

In B^ — > D*^D*~ decays the final state D* mesons may be in a state of 
S-, p- or d-wave relative orbital angular momentum. Since s- and d- waves 
are even under the CP transformation while the p-wave is odd, the CP- 
violation parameters in Eq. (1) are diluted. In order to determine the dilution, 
one needs to measure the CP-odd fraction. This can be accomplished with a 
time-integrated angular analysis. The BaBar collaboration has measured the 
polarization and CP asymmetries [8], and find the CP-odd contribution to 
be small, consistent with theoretical expectations [7]. The CP asymmetries 
are found to differ slightly from the expectation that neglects the contribution 
from the penguin diagram. 

In this Letter we report measurements of the branching fraction, the polariza- 
tion parameters and CP asymmetries in B^ D*^D*~ decays based on a 140 
fb~^ data sample, which corresponds to 152 million BB pairs. At the KEKB 
energy-asymmetric e^e~ (3.5 on 8.0 GeV) collider [9], the T(4S') is produced 
with a Lorentz boost of = 0.425 antiparallel to the positron beamline {z) . 
Since the P° and P° mesons are approximately at rest in the T(45') center-of- 
mass system (cms). At can be determined from the displacement in z between 
the fcp and /tag decay vertices: At ~ {zcp - Zug)/{P^c) = Az/{P^c). 

The Belle detector [10] is a large-sohd-angle spectrometer that includes a 
three-layer sihcon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters (ACC), time-of-flight 
(TOF) scintillation counters, and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron fiux-return located outside of the 
coil is instrumented with resistive plate chambers to detect mesons and to 
identify muons (KLM). 



2 Event Selection 



We reconstruct P° — > D*^D* decays in the following D* final states; (£)°7r+, 
D°7r-), (P)°7r+, P)-7rO) and {D+n^ ,Wt^-). For the decays we use 
K~Tr+, K^TT+TT^, K-n+n+n-, K+R-, K^tt+tt- and Kln+rr-n°. For the D+ 
decays we use D+ K^7r+, K^tt+tt'^, K-tt+tt+ and K-K+tt+. We 

allow all combinations of D decays except for cases where both D decays 
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include neutral kaons in the final state. 



Charged tracks from D meson decays are required to be consistent with orig- 
inating from the interaction point (IP). Charged kaons are separated from 
pions according to the likelihood ratio Pr/tt = + ^ij^)]^ where 

the likelihood function C is based on the combined information from the ACC, 
CDC dE/dx and TOF measurements. We require Pk/-k > 0.1 (0.2) for kaons in 
2-prong (4-prong) D meson decays. The kaon identification efficiency is 96%, 
and 13% of pions are misidentified as kaons. Candidate charged pions are re- 
quired to satisfy Pk/tt < 0.9, which provides a pion selection efficiency of 91% 
with a kaon misidentification probability of 3%. Neutral pions are formed from 
two photons with invariant masses above 119 MeV/c^ and below 146 MeV/c^. 
To reduce the background from low-energy photons, we require > 0.03 
GeV for each photon and p^^o > 0.1 GeV/c, where E^ and p^o are the photon 
energy and the 7r° momentum in the laboratory frame, respectively. Candi- 
date — > TT+TT" decays are reconstructed from oppositely charged track pairs 
that have invariant masses within 15 MeV/c^ of the nominal mass. A re- 
constructed is required to have a displaced vertex and a fiight direction 
consistent with that of a Kg originating from the IP. 



Candidate D mesons are reconstructed from the selected kaons and pions, 
and are required to have invariant masses within 6cr {3a) of the D meson 
mass for 2-prong (3- or 4-prong) decays, where a is the mass resolution that 
ranges from 5 to 10 MeV/c^. In this selection a is obtained by fitting the 
Monte Carlo (MC) simulated D meson mass. These (-D"*") candidates are 
then combined with tt"*" (tt^) to form candidates, where the IP and pion 
identification requirements are not used to select vr"*" candidates. The mass 
difference between D*~^ and (D^) is required to be within 3.00 (2.25) 
MeV/c^ of the nominal mass difference. We identify B meson decays using 
the energy difference AE = E^^ — E^^^^ and the beam-energy constrained 

mass Mbc = \/(-E'{5cam)^ ~ {pW^y^ where E^^^^ is the beam energy in the cms, 
and E^^ and p^"^*^ are the cms energy and momentum, respectively, of the re- 
constructed B candidate. The B meson signal region is defined as \ AE\ < 0.04 
GcV and Afbc within 3a of the B meson mass, where a is 3.5 MeV/c^. In order 
to suppress background from the e~^e~ — > uu, dd, ss, or cc continuum, we 
require H2/H0 < 0.4, where H2 (Hq) is the second (zeroth) Fox- Wolfram mo- 
ment [11]. After applying this requirement, we find that the contributions to 
the background from B~^B~, B^B^ and continuum are approximately equal. 
Figure 1 shows the Mbc and AE distributions for the B^ D*^D*^ candi- 
dates that are in the AE and Mbc signal regions, respectively. In the Mbc and 
AE signal regions there are 194 events. 
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Fig. 1. (Left) and (right) /S.E distributions for S° D*'^D*~ candidates within 
the AE (Mbc) signal region. Sohd curves show the fit to signal plus background 
distributions, and dashed curves show the background contributions that comprise 
B^B^, B^B^ and continuum events. 

3 Branching Fraction 

To determine the signal yield, we perform a two-dimensional maximum like- 
hhood fit to the M^^-AE distribution (5.2 GeV/c^ < Mbc < 5.3 GeV/c^ and 
\AE\ < 0.2 GcV). Wc use a Gaussian signal distribution plus the ARGUS 
background function [12] for the Mbc distribution. The signal shape param- 
eters are determined from MC. The background parameters are obtained si- 
multaneously in the fit to data. The AE distribution is modeled by a double 
Gaussian signal function plus a linear background function. We obtain shape 
parameters separately for candidates with and without D*'^ — > D'^n^ decays 
to account for small differences between the two cases. 

The fit yields 130 ± 13 signal events, where 20% include D*+ D+n° decays. 
To obtain the branching fraction B(-B° — > £)*+D*~), we use the reconstruc- 
tion efficiency and the known branching fraction for each subdecay mode. 
We obtain an effective efficiency of [1.06 ± 0.08] x 10^^ from the sum of the 
products of MC reconstruction efficiencies and branching fractions for each of 
the subdecays. Small corrections are applied to the reconstruction efficiencies 
for charged tracks, neutral pions and Kg mesons to account for differences 
between data and MG. 

We obtain 

B(B° D*+D*-) = [0.81 ± 0.08(stat) ± O.ll(syst)] x 10'^ (2) 

where the first error is statistical and the second is systematic. The result is 
consistent with the present world-average value [13]. 
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Fig. 2. Definition of the angles in the transversity basis. Angle 9ti and (ptr are defined 
in the D*^ rest frame (the lower plane), while 9i is defined in the D*~ rest frame 
(the upper plane). 

The dominant sources of the systematic error are uncertainties in the tracking 
efficiency (11%) and in the subdecay branching fractions (7%). Other sources 
are uncertainties in the fit parameters and methods (1%), in the reconstruction 
efficiencies of 7r° (2%) and (1%), particle identification (1%), polarization 
parameters (2%), the number of B mesons (1%), and MC statistics (1%), 
where each value in parentheses is the total contribution. 



4 Polarization 

The time-dependent CP analysis requires knowledge of the CP-odd fraction. 
To obtain the CP-odd fraction without bias, we must take into account the 
efficiency difference between the two CP-even components. Therefore, we per- 
form a time-integrated two-dimensional angular analysis to obtain the fraction 
of each polarization component. We use the transversity basis [14] where three 
angles 9i, 6tr and (ptr are defined in Fig. 2. The angle 6i is the angle between 
the momentum of the slow pion from the D*~ in the D*~ rest frame and the 
direction opposite to B momentum in the D*~ rest frame. The angle ^tr is 
the polar angle between the normal to the D*~ decay plane and the direction 
of fiight of the slow pion from the D*^ in the D*^ rest frame. The angle 0tr 
is the corresponding azimuthal angle, where 0tr = is the direction antipar- 
allel to the D*~ fiight direction. Integrating over time and the angle 0ti., the 
two-dimensional differential decay rate is 



1 (py 

F d cos 6*^6? cos 6*1 
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^ RiHi{cos6tr,cos6i), 

i=0,l± 



(3) 
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where i = 0, ||, or ± denotes longitTidinal. transverse parallel, or transverse 
perpendicular components, Ri is its fraction that satisfies 

Ro + Rii + i?± = 1, (4) 



and Hi is its angular distribution defined as 



i/o(cos 9tr, COS 9i) = 2 sin^ 9tr cos^ 9i, 

if||(cos6'tr, cos6'i) = sin^ ^^tr sin^ 6'i, (5) 
i7_L(cos 9tr, COS 9i)= 2 cos^ ^tr siu^ 9i. 

The fraction R± corresponds to the CP-odd fraction. 

Eq. (3) is affected by the detector efficiency, in particular due to the correla- 
tions between transversity angles and slow pion detection efficiencies. To take 
these effects into account, we replace i7i(cos^tr, cos^i) with distributions of 
reconstructed MC events T-[i{cos9t^.,cos9i), which are prepared separately for 
candidates with and without D*~^ — > D^tt^ decays as is done in the branch- 
ing fraction measurement. We also introduce effective polarization parameters 
R[ = €iRi/{eQRQ + €\\R\i + e±R±), where is a total reconstruction efficiency 
for each transversity amplitude. As a result, the signal probability density 
function (PDF) for the fit is defined as 

'^sig = R'iHi{cos 9tT, cos 9i). (6) 

i 

We determine the following likelihood value for each event: 

^ = /sig'^^sig + (1 - fsigj'Hbg, (7) 



where /gig is the signal probability calculated on an event-by-event basis as a 
function of AE and M^c. The background PDF Hhg is determined from the 
sideband region (5.20 GeV/c^ < Mbc < 5.26 GeV/c^, \AE\ < 0.2 GeV). A fit 
that maximizes the product of the likelihood values over all events yields 

R± = 0.19 ± 0.08(stat) ± O.Ol(syst), 

Ro - 0.57 ± 0.08(stat) ± 0.02(syst). (8) 

Figure 3 shows the angular distributions with the results of the fit. 

We study the uncertainties of the following items to determine the systematic 
errors: background shape parameters, angular resolutions, and slow pion de- 
tection efficiencies. Also included are a possible fit bias, MC histogram bin size 
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Fig. 3. Angular distributions of the S° D*^D*~ candidates in (left) cos ^tr and 
(right) cos 01 projections. In each figure, the dot-dashed, dotted and dashed lines 
correspond to longitudinal, transverse parallel and transverse perpendicular polar- 
ization components, respectively. The thin solid line is the background, and the 
thick solid line shows the sum of all contributions. The asymmetry in the cos^i 
distribution is due to the inefficiency for low momentum track reconstruction. 

dependence and misreconstruction effects. These systematic errors are much 
smaller than the statistical errors. 



5 CP Asymmetries 

We perform an unbinned maximum likelihood fit to the three dimensional Ai, 
cos^tr and cos^i distributions for S° — > candidates to measure the 

CP- violation parameters. 

The meson decay vertices are reconstructed using the D meson trajectory 
and an IP constraint. We do not use slow pions from decays. We require 

that at least one D meson has two or more daughter tracks with a sufficient 
number of the SVD hits to precisely measure the D meson trajectory. The /tag 
vertex determination is the same as for other CP- violation measurements [4]. 

The 6-fiavor of the accompanying B meson is identified from inclusive prop- 
erties of particles that are not associated with the reconstructed P" fcp 
decay [3]. We use two parameters, q and r, to represent the flavor tagging 
information. The first, q, is already defined in Eq. (1). The parameter r is an 
event-by-event, MC-determined flavor-tagging dilution factor that ranges from 
r = for no flavor discrimination to r = 1 for unambiguous flavor assignment. 
This assignment is used only to sort data into six r intervals. The wrong tag 
fractions for the six r intervals, (/ = 1,6), and differences between P° and 
P° decays, Aw;, are determined from the data; we use the same values that 
were used for the sin 20i measurement [4] . 
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The signal PDF is given by 



sig 



4tbo 



R'{Hi{cos6tr:,cos6i) 



i=0,||,J 



X 



1 -qAw + q{l - 2w) {A COS Am At + ^iS sin Am At) 



(9) 



where CP parity is +1 for i = and ||, and —1 for i —±. We assume univer- 
sal CP- violation parameters in Eq. (9), i.e. Sq = S\\ = S± and Aq = A\\ = A±. 
The distribution is convolved with the proper-time interval resolution function 
Rsig{At) [4], which takes into account the finite vertex resolution. 

We determine the following likehhood value for the j-th event: 



Pj = (1 - /oi) /[/signig(At')i?sig(At, - At') 

+ (1 - /sig)nkg(Ai')i?bkg(Ai, - At')] d(At') + /oiPoi(Aii), (10) 

where Poi(At) is a broad Gaussian function that represents an outlier com- 
ponent [3] with a small fraction fo\. The /gig calculation is explained in the 
previous section. The PDF for background events, Pbkg(At), is expressed as 
a sum of exponential and prompt components, and is convolved with i?bkg 
that is a sum of two Gaussians. All parameters in Pbkg(At) and -Rbkg arc de- 
termined by a fit to the At distribution of a background- enhanced control 
sample; i.e. events outside of the AE-M\^c signal region. We fix r^o and Ama 
to their world-average values [13]. The only free parameters in the final fit 
are »S and A, which are determined by maximizing the likelihood function 
L = Ylj Pj{Atj, cos 9trj, cos 9ij;S, A), where the product is over all events. 
The fit yields 



S = -0.75 ± 0.56(stat) ± 0.12(syst), 

A-- -0.26 ± 0.26(stat) ± 0.06(syst), (11) 

where the first errors are statistical and the second errors are systematic. These 
results are consistent with the SM expectations for small penguin contribu- 
tions. 

We define the raw asymmetry in each At bin by (A^g=+i — A^g=-i)/(-/Vq=_,_i -|- 
Ng=_i), where A^q=_|_i(_i) is the number of observed candidates with q — 
+1(— 1). Figure 4 shows the raw asymmetries in two regions of the flavor- 
tagging parameter r. While the numbers of events in the two regions are 
similar, the effective tagging efficiency is much larger and the background di- 
lution is smaller in the region 0.5 < r < 1.0. Note that these projections onto 
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Fig. 4. Raw — > D*^D*^ asymmetry in bins of At for (top) < r < 0.5 and 
0.5 < r < 1.0 (bottom). The solid curves show the result of the unbinned maxi- 
mum-likelihood fit. 

the At axis do not take into account event-by-event information (such as the 
signal fraction, the wrong tag fraction and the vertex resolution), which are 
used in the unbinned maximum- likelihood fit. 

The sources of the systematic errors include uncertainties in the vertex recon- 
struction (0.05 for S and 0.03 for A), in the flavor tagging (0.04 for S and 
0.02 for A), in the resolution function (0.05 for S and 0.01 for A), in the back- 
ground fractions (0.04 for S and 0.02 for A), in the tag-side interference [4] 
(0.01 for S and 0.03 for A), and in the polarization parameters (0.06 for S 
and 0.01 for A). Other contributions for S come from a possible fit bias (0.04) 
and from uncertainties in r^o and Arrid (0.02). We add each contribution in 
quadrature to obtain the total systematic uncertainty. 

We perform various cross checks. A fit to the same sample with A fixed at 
zero yields S = — 0.69 ±0.56(stat). We check with an ensemble of MC pseudo- 
experiments that the fit has no sizable bias and the expected statistical errors 
are consistent with the measurement. We also select the following decay modes 
that have similar properties to the 5° D*^D*~ decay: D*~D*~^, 
D-D*+, D*-Df, D'D+, and B+ D*°D*+, WD*+, D*^Df and WDf. 
Fits to the control samples yield S[B^ /:)(*)/}(*)] = _o.l2 ± 0.08, 

= +0.02 ± 0.05, S[B+ = -0.10 ± 0.07, and A[B+ 

= -0.001 ± 0.050, where errors are statistical only. All results are 
consistent with zero. We also measure the B meson lifetime using B^ 
D*+D*- candidates as well as the control samples. All results are consistent 
with the present world-average values. A fit to the At distribution of the B^ 
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£)*+£)*- without using polarization angle information yields <S = — 0.57±0.45, 
A = —0.29 ± 0.26; this result suggests that the CP-odd component is small, 
supporting our polarization measurement. 

Although the statistics are not sufficient to provide tight constraints, we also 
consider polarization-dependent values for S and A, which may arise from 
possible differences in the contributions of the penguin diagrams. We assume 
that the CP asymmetries for the CP-odd component are consistent with 
the SM expectations, and fix S± at the world-average value of sin20i [13] 
and A± at zero. A fit with this assumption yields S — —0.72 ± 0.50 and 
A = —0.42 ± 0.30 for the CP-even component, also consistent with the SM 
expectations. 



6 Conclusion 

In summary, we have performed measurements of the branching fraction, the 
polarization parameters and the CP-violation parameters for P° D*^D*~ 
decays. The results are 

B{B^ D*^D*-) = [0.81 ± 0.08(stat) ± O.ll(syst)] x 10-^ 
P^ = 0.19±0.08(stat) ±0.01(syst), 
Ro = 0.57 ± 0.08(stat) ± 0.02(syst), 
S = -0.75 ± 0.56(stat) ± 0.12(syst), 
^= -0.26 ± 0.26(stat) ± 0.06(syst). (12) 

The polarization parameters and CP- violation parameters are consistent with 
the SM expectations and theoretical predictions for small penguin contribu- 
tions [15]. 
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